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Big bang nucleosynthesis (BBN) is affected by the energy density of a primordial magnetic field 
(PMF). For an easy derivation of constraints on models for PMF generations, we assume a PMF 
with a power law (PL) distribution in wave number defined with a field strength, a PL index, and 
maximum and minimum scales at a generation epoch. We then show a relation between PL-PMF 
parameters and the scale invariant (SI) strength of PMF for the first time. We perform a BBN 
calculation including PMF effects, and show abundances as a function of baryon to photon ratio rj. 
The SI strength of the PMF is constrained from observational constraints on abundances of 4 He and 
D. The minimum abundance of 7 Li/H as a function of n slightly moves to a higher 7 Li/H value at a 
larger rj value when a PMF exists during BBN. We then discuss degeneracies between the PL-PMF 
parameters in the PMF effect. In addition, we assume a general case in which both the existence 
and the dissipation of PMF are possible. It is then found that an upper limit on the SI strength of 
the PMF can be derived from a constraint on 4 He abundance, and that a lower limit on the allowed 
7 Li abundance is significantly higher than those observed in metal-poor stars. 
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INTRODUCTION 



Element abundances in the early Universe are important observables to determine the features of the physical 
processes at the epoch of big bang nucleosynthesis (BBN). From a discovery of expansion of Universe [l| and the 
development of the cosmological theory [H 0, 0, 0; ; ^ was suggested that the early Universes was enough hot and 
dense enough to allow nuclear fusion to be operative, which is the original idea of BBN [a, LD, la l2| • ^ was assumed 
that radiative neutron capture reaction generates primordial elements in hot and dense conditions. However, neutrons 
and protons can be converted to each other via weak interactions. Hayashi studied this effect and proposed that BBN 
started from not a pure neutron gas, but a hybrid gas composed of neutrons and protons as well as radiation [loj . 
Hayashi studied routes of 4 He production and indicated that a relic abundance of 4 He depends on an initial abundance 
of the neutron [l(|. This suggestion and the consideration of relativistic quantum statistics refined the BBN model 
[ill ). The discovery of the cosmic microwave background (CMB) [l2j] indicated that the early Universes was hot and 
dense, and BBN theory was confirmed. Production of D, 3 He, 4 He, 7 Li [r|, all stable isotopes of Li, Be and B [l4j . 
and heavier elements up to 20 Ne 15| have been studied. Now the BBN theory includes many important physics, and 
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the numerical precision of the BBN calculations is rather high 

The baryon to photon ratio, 77, is one of the important parameters characterizing the BBN model. Recently, we 
have constrained rj from observations of CMB [HI, [3, E3, 0] and tested various models of BBN using the constrained 
value of 77 as an input. 

Effects of a primordial magnetic field (PMF) on BBN have been studied [49L [Hoi [5ll [H^ [53l [H^. [HBL [56l [57l [58L [59l| . 

The PMF directly affects rates of weak reactions [S^, [H(| , and also the cosmic expansion rate throu gh t he PMF en ergy 



density. Results of BBN could, therefore, be different from that for the case with no pmf [m, SHI Hi, 0, m 
The boost of cosmic expansion by the PMF energy density is the dominant effect of the PMF, and the change of weak 



reaction rates is subdominant [53l . 1 5 a . l56j ] . An updated result of light elements up to Li in the BBN model including 
a PMF has been derived based on a network calculation consistently taking account of important effects of PMF [H(| ■ 

A power law (PL) spectrum is one of the most familiar distributions for the energy densities of various physical phe- 
nomenon on cosmological scales, including the PMF(e.g. 5{J 61 , 62| and references therein). Thus many researchers 
have been studying the effect of the PMF with the PL spectrum on various physical processes in the early Universe 
5^ . 61, 62|,[63j] and have tried constraining parameters of such PMFs. 

In this article, providing the PL as the distribution function of the PMF (PL-PMF), we report effects of the PL- 
PMF on primordial abundances of light elements up to Li as a function of the baryon to photon ratio r\. The energy 
density of the PL-PMF is dependent on the field strength, the PL index, and the maximum and minimum scales of 
the PMFs at the generation epochs. Therefore, we also investigate degeneracies of these PL-PMF parameters in the 
effects of the PMF on light element abundances in detail for the first time. 

In Sec. II, we describe how to include effects of the PL-PMF on BBN in a code for the standard BBN (SBBN) model. 
In Sec. Ill, we release the numerical computational results of BBN with the PL-PMF, and we discuss comparisons 
between our results and the observational constraints on elemental abundances in Sec. IV. In Sec.V, we summarize 
our study and mention future works. 



MODEL 



In this section, we describe our code for the BBN calculation with recent updates, and suggest how to consider the 
effects of the PL-PMF on BBN. We only consider the effect of boosts of cosmic expansion by the PMF energy density 
because this effect is more important than other effects such as changes in the rate of weak 53, |55j, l56[ and nuclear 
[HoT ] reactions, and those in the distribution functions of electrons and positrons [13|, which were taken into account 
in Ref. [H(|. We modify the Kawano's BBN code 41, fiij, taking into consideration the PMF effects, and use the 
updated reaction rates with the recommendations of the JINA REACLIB Database VI. [HH as adopted in Ref. [60j . 



We use T n = 878.5(±0.7 s t ± 0.3 sy ) s as the neutron lifetime 661 ] . which is estimated with an improved measurement 
67j . This lifetime is somewhat shorter than the average neutron lifetime of 880.1 ±1.1 s calculated by the Particle 
Data Group [68} . It better satisfies the unitarity test of the Cabibbo-Kobayashi-Maskawa matrix [67} , and it slightly 
affects the primordial abundances (most significantly that of 4 He) [69(. Effects of the difference between this lifetime 
and a longer one 885.7 ± 0.8 s from the old recommendation by the Particle Data Group [7(| are seen in Figs. 1 
and 2 in Ref. [60]. The uncertainties in the calculated abundances originating from an uncertainty in the neutron 
lifetime are much smaller than both uncertainties in observationally induced abundances and changes in calculated 
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abundances caused by the magnetic field. We, therefore, do not take into account the uncertainty in the neutron 
lifetime. 



Constraints on primordial elemental abundances 

We shall describe constraints on the primordial elemental abundances that we adopt in this article. The primordial 
deuterium abundance is derived by observing the Lyman-a absorption systems of the quasi-stellar objects (QSOs). 
A recent measurement of the deuterium abundance in a damped Lyman alpha system of QSO SDSS J1419+0829 was 
performed more precisely than those of any other QSO absorption system [7l|. We then adopt two constraints: (1) 
a mean value of ten QSO absorption systems including J1419+0829, and (2) the abundance of J1419+0829 from the 
best measurement, as follows, 

2.40 x 10~ 5 < - < 2.88 x 10~ 5 
H 

(2cr range from mean value) , (1) 
D 
H 

(2<T range from best value). (2) 



2.43 x 10" 5 < - < 2.64 x 10 



The two groups IT10 [72j and AOS10 [73J have reported different results for the primordial abundance of helium, Y p , 
from surveys of metal-poor extragalactic HII regions. They obtained the following constraints, i.e., 

IT10 : 0.2463 < Y p < 0.2667 (2a), (3) 
AOS10 : 0.2345 < Y p < 0.2777 (2a). (4) 



Observations of the metal-poor halo stars (MPHSs) provide a constraint on the abundance of 7 Li 74j, i.e 



7 L : 

1.06 x 10~ 10 < — < 2.35 x 10~ 10 (2a). (5) 
H 



The Wilkinson Microwave Anisotropy Probe (WMAP) project, which observes the CMB across the full sky, con- 

0.00114 
0.00112 

0.309 



strains the baryon density parameter [48(, fl^h 2 = 0.02258lg 00112 (%&)■ This result corresponds to the baryon to 
photon number ratio, r^/lO" 10 = 6.225l^Jg (2a) 



Cosmic Expansion Rate and BBN 

In this subsection, at first, we introduce the cosmic expansion. Next we describe an effect of the cosmic expansion 
on BBN. Tthc Hubble parameter H for the homogeneous and isotropic flat-Universe is derived by the Friedmann 
equation and the conservation of energy momentum tensor, 

a\ 2 o 8ttG 

a) (6) 

f = - 3 "('+?). < 7 > 

where the dot expresses the time derivative as a = da/dt, G is the gravitational constant, c is the photon speed, p 
and p are the energy density and the pressure, respectively, in the Universe. The inverse of H corresponds to the age 
of the Universe, < agc ~ H . We assume that the main components of the Universe are photon, neutrino, baryon and 
cold dark matter (CDM), they are the ideal fluid, and the Universe expands adiabatically. 

We describe a rate of reaction between some particles as T. The time scale of the reaction is then When the 
reaction time scale becomes larger than the age of the Universe, i.e., i agc ~ H~ x , the reaction effectively stops in the 
Universe. Whether a reaction is operative or not is determined by the following condition, i.e., 

r _1 < H , r > H : the reaction continues, (8) 
r _1 > H' 1 , r < H : the reaction stops. (9) 

Since the final nuclear abundances of BBN depend on when the respective reactions become ineffective, the energy 
densities of the Universe influence the abundances of the elements produced in BBN. 
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Energy density of power law PMF 



We will explain how the PL-PMF energy density affects the cosmic expansion and BBN. In the early Universe 
including the BBN epoch, we can assume that statistically average motions of fluids are negligibly small. In this case, 
from Maxwell's equations and Ohm's law, we derive the induction equation [75| 

f = < 10 > 

Here £ is the magnetic diffusivity, which is defined by 

i i 

C C U) e , . 

where a is the electric conductivity and w e is the electrical resistance [TH as follow, 

1 m e ca T m e 

uj e = — = -criyCTT = 5— (12) 

a n e e z r\e* 

Here m e is the mass of the electron, e is the electron charge, <tt is the Thomson scattering cross section, n 7 is 
the number density of the photon, n e is the number density of the electron, and 77 is the baryon to photon ratio. 
Equation (fTU|) describes the dissipation of the magnetic field, and it shows that the magnetic field dissipates rapidly 



with time on the scale L < y/t^jT)^ [59|, l75[ , where T is the temperature, and t age is the age of the Universe which 



is defined by i agc ~ 1.32 x (lMeV/T) s at the relativistic degrees of freedom, = 3.36. Therefore, the magnetic 
field is hard to be generated on the scale L <C \/t agc (T)( = Lpj(T). On the other hand, the magnetic field on a scale 
much larger than Lpi(T) is hard to dissipate in the Universe by time t agc , and the magnetic field is "frozen-in" the 



dominant fluids [75[. In this article, we are interested in magnetic fields which survive until the end of BBN. In the 
epoch during which the electron positron annihilation has been completed (T ~ O.OlMeV), the cutoff scale for PMF 
generations, Lfi(T) for 10~ ln < 77 < 10~ 9 , is 

2. 7 3xlQ 5 cmx( - y cV ) i <L FI (r) 

< 8.63 x 10W (^^)\ (13) 

the comoving scale L F1 = L F1 {T = 0.01MeV)/a(T = O.OlMeV) is 

1.16 x 10 13 cm < L Fl < 3.67 x 10 13 cm (14) 

and the comoving wave number kpj — 2ir / L F1 is 

1.72 x lCT^cm" 1 < k FI < 5.41 x lCT^cm" 1 . (15) 

Here we use a(T = O.OlMeV) = T CM B/0.01MeV - 2.35 x 10~ 8 . Thus we will report results of the effects on BBN 
from the PMF with a comoving wave number less than fcpi ~ 10~ 13 cm -1 . 

The PL spectrum is one of the most familiar forms for distributions of the energy densities including that of the 
PMF [HI, |6l|, |62| |63j|. It appears in various physical phenomena operating on the cosmological scale. Many researchers 



have been studying the effect of the PMF with the PL spectrum on various physical processes in the early Universe 
[Hil WH [H^, H^], and have tried to constrain parameters of such PMFs. We then assume that the PMFs with the PL 
spectrum are statistically homogeneous, isotropic, and random, as assumed in previous studies, and we research their 
effects on BBN. 

After a magnetic field is generated, the total energy density originating from that of the magnetic field is conserved, 
even if part of the magnetic field energy is transferred to other physical degrees of freedom (e.g. kinetic or/and thermal 
energies of dominant fluids or the gravitational wave) in the adiabatic Universe. The cosmic expansion depends on 
the total energy density in the Universe. Therefore, when we consider the effect of the PMF energy density on the 
cosmic expansion, it is reasonable to define the total energy density originating from the PMF at its generational 
epoch. We shall discuss cases of finite energy transfer from the PMF to other physical degrees of freedom in Sec. „ 

Since the strength of the frozen-in PMF is proportional to an inverse square of the scale factor, its energy density 
and pressure depend on the inverse fourth power of the scale factor as follows [61], [62 , 



3 B 2 (a ) 

PPMF = ^-PPMF OC -. . (16) 
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We can neglect the pressure and energy density of matter in the radiation dominated era. Therefore, the total energy 
density and pressure in the early Universe are 

P ~ Pr + PPMFj (17) 

P~Pr + Ppmf- (18) 
Next, we shall derive the energy density of the PL-PMF. A two-point correlation function of the strength of the 
PMF [Hl|, [H[63, [Zl [Zllzl is given by 

{B\k)B *{U)) = (27r) 3 P [PMF] (k)Pi(k)5(k-k') . (19) 

Here 

P^) = 5)- k ^L (20 ) 

and 

P [PM F](k) = Ak n *. (21) 

The convention for the Fourier transform is 

f(k) = J cxp(ik- x)F{x)d 3 x. (22) 

In this article we need the energy density of the PMF in the Universe, ppmf = (B 2 )/(Stt). From Eq. (|19|) . we obtain 

(B l {k)BC{ti)) = 2(2n) 3 P [PMF] (k)S(k-k') 

= 2(2ir) 3 Ak nB 5(k- k'). (23) 

We derive the coefficient from the variance of the magnetic fields A in real space. First we define 

(B^B^l^Bl (24) 

Here A is a comoving radius for a Gaussian sphere, and B\ is a comoving PMF strength which is scaled to the 
present-day value on some length scale A. From Eqs. (jT9|) - (|2Tj) and ([24]) . 

(B i (x)B i (x))\ x = Bl 

= - \ e f d 3 k f d 3 k' exr>(-ix- k + ix- k') 
(27r) 6 7 J 

x(B\k)B:(k'))x\Wi(k)\, (25) 
where W\(k) is a window function here assumed to be a Gaussian filter, W\{k) = exp(— X 2 k 2 /2). We then obtain the 



coefficient A as follows 



A = B 



(2tt 



,2 



2 



J dkk nB+2 exp(-A 2 fc 2 ) 



where A = 27r/fcpMF]- Therefore 



B2 (2n) 2 A"b+3 

^ 2 r ^ " b + 3 ) 

[PMF] V 2 ; 



CO 7r ^2 R 2\n B +3 



From Eqs ([191) " (ESJ) and ([27]), the energy density of the PMF is 

?2\ o /./cr maxl J7 7.3 



2 \n B +3 



2 fWl dfc fc 3 {2n) 2 B 2 \ nB + 3 



Stt J k[mia] k 2tt 2 2T (2iLfc3) 
J ^— [(Afc [max] )" B + 3 - (Afc [min] )" B + 3 ] 



87rr(2^ 



(28) 
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where fc[ ma x] and fc[ m j n ] are the maximum and minimum wave numbers, respectively, which depend on generation 
models. 

The first purpose of this study is to investigate the effects of the PMF energy density on BBN and to constrain the 
PMF energy density in the BBN epoch. The second purpose is to discuss of the degeneracy between the PL-PMF 
parameters and distribution models of the PMF. To effectively proceed with these studies and discussions, we redefine 
the scale invariant (SI) field strength of the PMF from Eq. (|2"51) , which is constrained from observational constraints 
on elemental abundances, i.e., 



R X _ R X [(Afc [m a x])" B+3 

^SI = \ \ 



(Afc [n 



\riB+3 



r(a*t£) 



(29) 



and 



B X (riB, fe[max], fymin]) = B 



fe „B+3 _ k n B + a ) A „ B - 
max mm 



<„n B +3 



(30) 

where X is a species of light element produced in BBN and -Bgj represents the corresponding value of the SI PMF 
strength. Since B^ is proportional to -^/ppmf and does not depend on other PMF parameters, it is useful to directly 
understand effects of the PL-PMF amplitude on BBN. Furthermore it is relatively easy to discuss the degeneracy of 
the PMF parameters if we use this formalization. 



RESULTS 



In this section, we show calculated results of the BBN with the PL-PMF using our code. We also show how the 
strength of the PL-PMF is constrained from limits on abundances of light elements up to Li. 



BBN with PL-PMF 

Figure Q] shows abundances of 4 He (Y p : mass fraction), D/H, 3 Hc/H, 7 Li/H and 6 Li/H as a function of r\ for cases of 
BBN with (solid lines) and without (dashed lines) the effects of the PL-PMF. The parameters of the PL-PMF are set 
to be B\ m P c = 2.2 fiG, iib = —2.99, fcr max ] = 10~ 13 cm -1 , and fcr min i = cm -1 . This parameter set corresponds to a 
SI PMF strength of i?gj = 2.54 /xG. The vertical painted band shows the limit on the baryon to photon ratio derived 
from WMAP 7yr data 48]. The horizontal painted bands and lines with downward arrows indicate observational 
constraints and upper limits, respectively, on elemental abundances. 

Effects of the PMF on abundances of respective nuclides have been analyzed in detail [&cj . The effects are described 
below to explain the abundances predicted as a function of r\ in Fig. [T] 

When the PL-PMF exists, abundances of the neutron are larger than those of the BBN without the PL-PMF, 
i.e., the SBBN. The reason is as follows: Considering the energy density of the PL-PMF, the cosmic expansion rate 
H is larger than that of SBBN. Therefore, the weak reaction freeze-out occurs earlier [5|| [cf. Eqs. <JSj> and (j9|)], 
and abundances of the neutron at the freeze-out increase. Furthermore, the faster expansion leads to a shorter time 
interval between the freeze-out and 4 He production, so abundances of the neutron, which is spontaneously decaying 
after the freeze-out, are higher at the epoch of 4 He production. Since 4 He nuclei are formed from almost all neutrons, 
the energy density of the PL-PMF increases the abundance of 4 He [Fig. [Ha)]. 

Deuterium ( 2 H, D) is produced and destroyed by 

p + n^ 2 H + 7 , (31) 
2 H+ 2 H^ 3 H+p, (32) 
2 H+ 2 H^ 3 Hc + n. (33) 

The high neutron abundance in the PL-PMF case leads to an increase in D abundance because of an enhanced 
production rate through the reaction (l3~Tj) [Fig. QJb)]. 

3 H and 3 He are produced via the reactions (f32|) and (l33l) . and destroyed by 

3 H + 2 H -> 4 He + n, (34) 
3 He + n^ 3 H+p. (35) 
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The enhanced abundances of the neutron and D result in higher 3 H and slightly higher 3 He abundances than those 
of SBBN. 3 H nuclei eventually j3~ decay to 3 He. The final primordial abundance of 3 He is then the sum of those 
of 3 H and 3 He in the BBN epoch. Primordial 3 He is mainly produced as 3 He (in SBBN more than 99.5 % of the 
final abundances of 3 He result from those of 3 He for ij — 10 -9 -10~ 10 ). The slight increase in 3 He at the BBN epoch 
enhances the final 3 He abundance slightly [Fig. Hfb)] . 



7 Li is mainly produced and destroyed by 



4 He + 3 H — > 7 Li + 7, (36) 
7 U+p^ 4 Hc+ 4 He. (37) 



A combination of a significant increase in 3 H abundance and a small fractional increase in 4 He abundance, both by 
the PL-PMF, results in a significant increase in the 7 Li production rate [Eq. (13T51) ]. while a small fractional decrease 
in H abundance (its mass fraction is X p 1 — Y p ) results in a slight decrease in the 7 Li destruction rate [Eq. (|37|) ]. 
As a result, 7 Li abundance is larger than that in SBBN. 
7 Be is mainly produced and destroyed by 

4 Hc + 3 He 7 Be + 7, (38) 
7 Be + rw 7 Li+p. (39) 



A significant increase in neutron abundance by the PL-PMF results in an increase in the 7 Be destruction rate [Eq. (|39|) ] , 
while a combination of a slight increase in 3 He abundance and a small fractional increase in 4 He abundance results 
in a slight increase in the 7 Be production rate [Eq. (f3"5|) ]. As a result, 7 Be abundance is smaller than that in SBBN. 

7 Be nuclei eventually capture electrons and are converted to 7 Li nuclei at the epoch of cosmological recombination. 
The primordial 7 Li abundance is, therefore, contributed from both 7 Li and 7 Be produced at BBN. If the baryon to 
photon ratio ij is lower than some critical value ?7 cr it, abundance of 7 Li is larger than that of 7 Be. On the other hand, 
if 7] is larger than ?7 cl it, abundance of 7 Be is larger than that of 7 Li. Therefore, the equality, i.e., 7 Li/H= 7 Bc/H, 
is satisfied at r\ — r/dit- Since the PMF increases 7 Li abundance, while decreases 7 Be abundance, the final 7 Li 
abundances are higher in rj < ry cr it and lower in r\ > r? cr ;t than those in SBBN. It should be noted that the 77 cr ;t value 
and the minimum value of the final 7 Li abundance, which is located at a valley of the curve for 7 Li/H, are moved 
from those of SBBN by the PMF effects: The values of SBBN are ?7 crit = 3.00 x 10~ 10 and 7 Li/H min = 1.18 x 10~ 10 
at 77 = 2.50 x lO" 10 , while those of BBN with the PL-PMF are 77 crit = 3.71 x lO" 10 and 7 Li/H min = 1.41 x lO" 10 at 
X) = 3.10 x 10- 10 [Fig.QJc)]. 



Li is produced and destroyed by 



4 He + 2 H — > 6 Li + 7, (40) 
6 Li+p^ 3 He+ 4 He. (41) 



The relatively large enhancement of D abundance and a small fractional increase in 4 He result in an efficient production 
of 6 Li [Fig. OK)]. 

The 4 He abundance increases as a function of 77, and it increases as a function of the PMF amplitude. The limits 
on 4 He abundance, therefore, prefer an 77 value smaller than in the SBBN when the effect of the PMF is considered. 
The D abundance is smaller for larger 77 values, while it is larger for larger amplitudes of the PMF. The limits on D 
abundance then prefer an 77 value larger than in the SBBN. The 7 Li abundance is larger for larger 77 values, while it 
is smaller for larger amplitudes of the PMF. The limit on 7 Li abundance, therefore, prefers an 77 value larger than in 
the SBBN when the effect of the PMF is considered (see Sec.„ also). 



Constraint on the PL-PMF from light element abundances 

To research not only qualitative results but also the quantitative relationship between Bgj and 77, we show how 
the PL-PMF strength Bfj can be constrained via observational constraints on the abundance of light elements 
[Eqs. (LT]) (J2]), d3j), d4j) and ([5])] and the 77 value. Figure [2] shows calculated constraints of B^s as a function of r\. 
The vertical painted band is the limit on the baryon to photon ratio from the WMAP 7yr data (48fl . We use abun- 
dances of only 4 He, D, and 7 Li to constrain the B^ Y value since limits on Bgj from abundances of 3 He and 6 Li are less 
severe than those of 4 He, D, and 7 Li [Figs. [TJb) and (d)]. Four kinds of lines in Fig. [2ja) correspond to upper limits, 
i.e., Bf v from Y p < 0.2777 (AOS10; black-thin-dotted-dashed line), Y p < 0.2667 (IT10; green-thick-dotted-dashed 
line), D/H < 2.88 xl0~ 5 (mean; blue-thin line), D/H < 2.64 xl0~ 5 (best; red-thick line), and 7 Li/H < 2.35xl0" 10 
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(purple-dotted line) , respectively. On the other hand, requiring a consistency between 7 Li abundances of theory and 
observations leads to a lower limit on the PL-PMF as a function of 77 [the right purple-dotted line in Fig.^a)] because 
the abundance of 7 Li decreases when the energy density of the PL-PMF increases (Sec,). Then, the allowed region 
from the 7 Li limit on the B S Y vs. rj plane is between the two purple-dotted curves in Fig. [Ha). 

The limits on PMF strengths, i.e, Bg\, at the r\ value of a 2 a lower limit by the WMAP 7yr data [48j], i.e, 77=5.916 
xlO~ 10 , are 

b\\ < 2.00 fiG for Y p < 0.2667 (IT10), (42) 
B^ < 2.56 fiG for Y p < 0.2777 (AOS10), (43) 
Bsi < 0.848 fiG for D/H < 2.88 x 10" 5 (mean), (44) 
B 7 S Y > 5.05 fiG for 7 Li/H < 2.35 x 10~ 10 , (45) 
at t) = 5.916 x 10~ 10 . 

Since the abundance of D/H is already more than D/H = 2.64 x 10~ 5 (best) when no PMF is assumed, we cannot 
put a constraint (Sg) using D/H < 2.64 x 10~ 5 at 77 = 5.916 x 10~ 10 , 

When we adopt the 77 values of the best-fit and 2 a upper limit by the WMAP 7yr data [48], constraints from 
7 Li abundance are weaker, or B S Y values are larger. We then show four constraints (-Bgj), except for that of 7 Li as 
follows: 

Yp< 0.2667 (IT10), (46) 
Yp< 0.2777 (AOS10), (47) 
D/H < 2.88 x 10 -5 (mean), (48) 
■ D/H < 2.64 x 10~ 5 (best), (49) 
at r? = 6.225 x 10~ 10 , 



B SI 


< 


1.97 fiG 


-°SI 


< 


2.54 fiG 


Bsi 


< 


1.51 /jG 


B SI 


< 


0.780 fxG 



and 





< 


1.95 fiG 


for 


Y p < 0.2667 (IT10), 


v 

D P 

B si 


< 


2.51 ^G 


for 


y p < 0.2777 (AOS10), 


B si 


< 


1.97 jiG 


for 


D/H < 2.88 x 10~ 5 (mean), 


B SI 


< 


1.45 fj,G 


for 


D/H < 2.64 x 10~ 5 (best), 



(50) 
(51) 
(52) 
(53) 



at r; = 6.539 x 10~ 10 . 



This result is very similar to the latest result [60j of the BBN calculation which consistently includes effects of the 
PMF through changes in the distribution functions of electrons and positrons as well as the cosmic expansion rate. 
This is because the effect of the PMF on BBN through changes in distribution functions is smaller than that through 
the cosmic expansion rate (56| . 

Adopting our results [ (|42|) -(|53 |) and Fig. [2] and the WMAP-77 value [48(, we show combined constraints on the 
PL-PMF [Fig. [UJb)]. If we adopt limits on abundances of F p (IT10) and D/H (mean), we obtain a combined upper 
limit: 

Bs! +B < 1-95 fiG 

forFp < 0.2667 (IT10) 

and D/H < 2.88 x 10" 5 (mean). (54) 

The maximum value corresponds to the limit at 77 = 6.523 x 10 -10 (see the intersection of the thin solid curve and 
the thick dot-dashed curve in Fig. [5]). If we adopt Y p (IT10) and D/H (best), the following upper limit is derived: 

Bli +D < 1-45 fiG 

forYp < 0.2667 (IT10) 

and D/H < 2.64 x 10^ 5 (best). (55) 

The maximum value corresponds to the limit at 7/ = 6.539 x 10 -10 (see the intersection of the thick solid curve and 
the vertical painted band in Fig. [2]). 
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DISCUSSIONS 



In the above sections, we explained the effect of PL-PMF on BBN and derived constraints on PMF energy 
densities. In this section, considering these results, we will discuss degeneracy between the PL-PMF parameters, 
(B\, tt-b, fc[max]i fymin]) i n effect s on BBN, how such constrained PL-PMF parameters are applied for constraint on 
PMF generation models, and cases of an energy transferring from PMF to others physical degrees of freedom. 



Degeneracy between PL-PMF parameters 

The PL-PMF strength Bf depends on fc[ max ], fc[ m in] an d n B [Eq. (J5Q]> ] ■ In particular, if fc[ m i n ]/fc[ max ] is close to 1 
or/and ns is close to -3 [80j], the influence of fc[ min ] on cannot be ignored and we should treat the fcr m ; n i value 
carefully in constraining PMF parameters and generation models (Fig. [3]), although fcr m i n i does not affect the energy 
density of the PL-PMF otherwise. If fc[ m i„]/fc[ max ] is close to 1, such a PMF is generated on narrow scale ranges. In 
this case, the distribution of the PMF may be better described by the Gaussian, Delta or log-normal |81 J functions 



than the PL. When tib is close to -3, the PMF is expected to be generated in the epoch of inflation [6f|, |62|, |79|. In 
this case, from Eq. (|30l) and Fig. [3] (b), the term of k"^^ cannot be neglected even if fc m i n /fc max is relatively small. 
Thus we need to be careful in constraining the PMF of inflationary origin. 

We have discussed the specific cases of the PL-PMF. Next, we discuss the opposite cases in which ranges of wave 
number for the generated PL-PMF spectrum are wide enough so that fc[ m i n ] can be neglected (fcr m j n ] /fcr roax i -C 1). 
Figure 0] shows constraints on the field strength B\ as a function of ub [cf. Eqs. (f30|) and (|54|) ] for fixed values of 
A = I Mpc and fc[ m i n ] = 0. It is clear that the limits on this plane of B\ vs. iib from light element abundances, i.e., 
their negative slopes, show a negative correlation of the parameters [Eq. ([30p ] in a parameter region for the fixed SI 
strength of the PMF. When generated PL-PMFs have spectra with smaller k\ max i, constraints on the (B\, ns) plane 
are more stringent. In general, considering the cosmic expansion after the inflation, the smaller the generation scale 
of the magnetic field, the earlier the magnetic field is generated. This feature is useful in constraining models for 
PMF generation in the post-inflation epoch. 

Since cosmological observations have constrained strengths of the PL-PMF B\ at A = I Mpc that are less than 



orders of I nG [62|, |73, |82|, |83j, |8J] , we shall discuss the degeneracy between /q max ] and ub in the PMF effect on BBN 
for 0.1 nG < B\ < 10 nG. Figure [5] shows constraints on the maximum wave number fc[ max ] as a function of tib [cf. 
Eqs. ((30)) and (|54)) ] for fixed values of A = I Mpc and /c[ m i n ] = 0. Results for various B\ values of 10 nG (dashed line), 
I nG (solid line) and 0.1 nG (dot-dashed line) are shown. It is obvious that the limits on the plane of fc[ max ] and n-g, 
show a negative correlation of the parameters [Eq. (|30|) ] in a parameter region for the fixed SI strength of the PMF. 



The cosmological observations (e.g. the CMB and the large scale structures: see Refs. [62J, 113, |82j, l£3j, |84{ for details) 
have constrained the PMF with strength of the order of I nG at I Mpc. As mentioned in Sec. „ we are interested in 
the magnetic field on scales much larger than the cutoff scale for the PMF which is generated before the end of the 
BBN epoch, fcpi ~ I0~ 13 cm _1 . In the case where ub is less than -2.6 and B\ = 10 nG, the fc[ max ] value is allowed to 
be larger than 10~ 13 cm _1 (Fig. [5]). On the other hand, if tib is larger than -2.6, the allowed value of fcr max ] becomes 
smaller. However the upper limits on fc[ max i for very large tib are smaller than the cosmic expansion rate, i.e., H (the 
lower limits on fcZ^i are larger than the Horizon scale) before BBN, and such a PMF can only be generated before 
the end of inflation. Therefore, using this feature, the PL index of the PMF generated in the post-inflation epoch can 
be strongly constrained. 



Energy transfer from PL-PMF 



From Fig. [21 it is confirmed that this BBN model including the effect of the PMF energy density cannot be a solution 
to the discrepancy between 7 Li abundances produced in SBBN and those deduced from observations of MPHSs. We 
consider a case where part of the PMF energy (pd) is dissipated and transferred to the radiation energies. The final 
value of the photon energy density measured at an epoch of scale factor a is then p 1 {a) = p^{a) + pd(a), where p^(a) 
is a primary photon energy density which does not include that transferred from the PMF by dissipation. Since 
Pd/p-y increases as a function of time from the time of the PMF generation, the following equation always holds: 
Pd (°bbn)/ /?<y(obbn) < Pd(oo)/ p-/(ao). This means that the baryon to photon ratio at the BBN epoch is expected to 
be more than the present value, and the rj value in the BBN epoch constrained from the WMAP 7yr data would be 
larger than the case of no energy transfer. 
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Here, we conservatively consider the most general case in which the history of generation and dissipation of the 
PMF in the Universe is completely unrestricted. The amplitude of the PMF in the BBN epoch and the largeness 
of the PMF energy dissipation during a time from BBN to the cosmological recombination should then be regarded 
as two independent parameters. The effect of the PMF on the cosmic expansion rate increases the abundances of 
4 He, D and 3 He, and decreases that of 7 Li for r\ values around the WMAP 7yr estimation for the SBBN (?7wmap)- 
The effect of dissipation, on the other hand, increases abundances of 4 He and 7 Li and decreases those of D and 3 He 
(Figs. [T]and [2]). Effects of the PMF and dissipation on the abundances of D, 3 He and 7 Li work in opposite directions, 
while they both increase 4 He abundance. One can, therefore, understand that an upper limit on the strength of the 
PMF from He, i.e., B s f , can be derived by assuming that a PMF dissipation does not occur. In this regard, the 
case of Fig. [T] corresponds to the maximum amplitude of the PMF allowed by the constraint from 4 He abundance for 
V > Iwmap- The minimum 7 Li abundance in the BBN model with a free parameter of B^ and 77 > ?7wmap is given 
by the cross region of the solid line for 7 Li/H and the vertical band (Fig. [T]), which is higher than the abundances in 
MPHSs. We then conclude that the abundances of 4 He and 7 Li cannot be simultaneously consistent with observations 
[Eqs. ([3]) and ((SJ)], even in the case where both the energy density and the dissipation of the PMF are taken into 
account. 



SUMMARY 



We showed quantitative and qualitative effects of a PMF on each abundance of a light element using a BBN model 
taking account of the PMF with two parameters of the SI strengths of PMFs and the baryon to photon ratio 77 (Figs. 
CD and©. 

We assumed a PMF characterized by the power spectrum of a PL in wave number. Parameters of this PL-PMF 
are a field amplitude B\, a PL index ne, and maximum and minimum scales at a generation epoch fcr max i and fer m i n ], 
respectively. We then showed a relation between PL-PMF parameters and the SI strength of the PMF which is the 
parameter constrained from observations of light element abundances (Sec. p). 

We then calculated BBN taking into account recent updates on reaction rates and the effect of the PMF, and 
showed results of the abundances in the example case of (Bi Mpc, fc[max]i ^[min]) = (2-2 /iG, —2.99, 10~ 13 cm" , 
cm -1 ) as a function of 77. It was observed that the critical 77 value, i.e., ?7 cr it, at which 7 Li and 7 Be are produced in 
the same abundance at the BBN epoch, moves by PMF effects. The value 77 cr jt = 3.00 x 10~ 10 for SBBN is enhanced 
up to 7i cr ;t = 3.71 x 10~ 10 for the example case. In addition, the minimum of the final 7 Li abundance, located at 
the valley of the 7 Li/H vs. r) curve moves from 7 Li/H rain = 1.18 x 10~ 10 at 77 = 2.50 x 10~ 10 for the SBBN to 
7 Li/H min = 1.41 x 10~ 10 at 77 = 3.10 x lO" 10 for BBN with the PL-PMF (Sec. 

We investigated constraints on parameters of the PL-PMF from the abundances of light elements up to Li produced 
in BBN. We showed respective and combined constraints on the parameters. As a result, we obtained upper limits 
on the SI strength of the PMF, i.e., B§f +D < 1.45 — 1.95/1 G from abundances of 4 He and D for the 2 a region of 77 
values by the WMAP 7yr data (Sec. J. 

In Sec. IV, we discussed the degeneracy of the PL-PMF parameters in effects on BBN and showed some possibility 
of constraining models of PMF generations by combining constraints on the SI field strength from light element 
abundances to the relation between PL-PMF parameters and the SI field strength (Sec. ^. 

In addition, we consider a general case in which the existence and energy dissipation of the PMF are allowed. Based 
on our result of the BBN calculation (Figs. [I] and [2]), it was found that an upper limit on the strength of the PMF can 
be derived from a constraint on 4 He abundance. A lower limit on 7 Li abundance is also derived, and it is significantly 
higher than those observed in metal-poor stars. We then conclude that it is impossible to solve the 7 Li problem by the 
PMF energy density, even if we consider that part of the PMF energy is dissipated and transferred to the radiation 
energies (Sec.,). 

In the near future, cosmological observations of the CMB and large scale structures will make it possible to constrain 
the PL-PMF parameters more tightly. If we combine such results and our constraint on the parameters through the 
study of BBN, we will not only limit the PL-PMF parameters but also the generation models of the PMF more 
precisely. 

We are grateful to Professor G. J. Mathews and Professor T. Kajino for valuable discussions. This work was partly 
supported by Grant-in- Aid for JSPS Fellows No.21.6817 (M.K.). 



11 



BBN without PL-PMF 

BBN with PL-PMF, 1 

(B x , « B , k [m . dx] , k [aiD] ) = (2.2 nG, -2.99, 10" 13 cm" 1 , cm" 1 ) J ~ 

5.916 <r|X10 10 < 6.539 H 




FIG. 1. Abundances of Y p (mass fraction of 4 He), D/H, 3 He/H, 7 Li/H and 6 Li/H as a function of rj in BBN models with and 
without the effect of the PL-PMF for the fixed value of A = 1 Mpc. The vertical painted band is the limit on the baryon 
to photon ratio from WMAP 7yr data [48l |. The white dashed line in the limited range of r\ indicates the best value. The 
horizontal painted bands and lines with downward arrows indicate observational constraints and upper limits, respectively, on 
elemental abundances. The constraints on 3 He and 6 Li abundances are from Ref. [85[ and Refs. [8a. l87j . respectively. 
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